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Abstract.  A  semi-classical  model  of  the  reflectance  of  a  complex  material  is  modeled.  Model  selection  techniques  using 
this  reflectance  model  are  developed  to  detect  molecular  vibration  modes  measured  using  infrared  spectroscopy.  An  inverse 
problem  methodology  is  then  used  to  ascertain  the  underlying  vibrational  modes  associated  with  varying  concentrations  of 
different  chemical  species.  These  two  techniques  are  developed  to  characterize  the  extent  of  material  degradation  through  the 
detection  of  spectra  associated  with  the  chemical  processes  underlying  the  degradation  process.  These  methods  are  validated 
using  data  from  samples  of  ceramic  matrix  composites.  The  samples  considered  were  exposed  to  varying  amounts  of  thermal 
treatment  to  ascertain  the  effectiveness  of  these  methods  to  detect  degradation. 

Keywords:  Reflectance  Spectroscopy,  Optical  Constants,  Fourier  Infrared  Spectroscopy,  Data  Analysis 


INTRODUCTION 

Evaluation  of  ceramic  matrix  composites  is  the  main  motivation  for  this  study,  however,  we  hope  that  these  methods 
or  methods  like  these  could  eventually  be  applied  to  other  materials.  Degradation  due  to  oxidation  takes  place  in 
ceramic  matrix  composite  before  an  actual  damage  occurs.  This  degradation  causes  these  materials  to  become  brittle 
which  leads  to  catastrophic  failure. 

Ceramic  matrix  composite  have  a  multiscale  geometry  which  causes  exact  modeling  of  the  components  to  be 
nonviable.  The  mesoscale  structure  includes  tows  which  are  woven  into  cross  plies  through  the  matrix;  this  is  on  the 
1  mm  scale.  The  tows  are  composed  of  fiber  bundles  which  are  on  the  scale  of  100  jum.  An  approximate  description 
is  depicted  in  Figure  1.  The  degradation  that  we  are  concerned  with  occurs  at  the  scale  of  the  individual  fibers  (~10 
nm),  see  Figure  2.  For  nondestructive  evaluation  applications,  the  exact  geometry  is  not  available  a  priori,  would  be 
too  computationally  intensive  to  include  in  models,  and  is  not  of  primary  interest.  Models  of  very  different  forms  can 
be  deduced  from  various  and  assorted  assumptions  though  none  of  these  can  expected  to  be  exact  [1].  In  our  analysis, 
we  used  model  selection  results  to  make  data  driven  choices  between  two  models  that  we  used  to  analyze  the  data. 

A  further  complication  in  this  problem  is  related  to  nature  of  the  degradation  that  we  seek  to  characterize.  At  high 
temperatures,  air  causes  the  formation  of  oxides  in  ceramic  matrix  composites.  This  oxidation  at  small  levels  (before 
it  is  observable  by  most  other  methods)  causes  the  material  to  be  brittle  [2].  The  degradation  can  be  observed  by 
transmission  electron  microscopy  and  electron  energy  loss  spectroscopy  [3],  These  methods  are  not  capable  of  in  field 
inspection  so  we  seek  a  method  that  could  be  performed  in  field.  Fourier  transform  infrared  (FTIR)  devices  have  been 
developed  with  in  field  inspection  capabilities  [4], 

Spectroscopy  is  a  very  common  technique  in  many  fields.  It  is  a  ubiquitous  tool  in  control  processes  in  the 
food  production,  medicine  and  petroleum  production,  to  name  a  few  [4].  It  is  also  used  in  detection  in  chemistry 
and  astronomy  [5],  Many  of  the  fields  consider  applications  where  Beer’s  law  applies  which  allows  the  data  to  be 
transformed  to  a  linear  model.  In  these  cases,  all  the  techniques  for  linear  models  may  be  applied  such  as  PCA  and 
PLS  [6,  7].  These  techniques  can  be  used  for  data  from  nonlinear  models  however  the  results  will  not  correspond  to  a 
physically  meaningful  quantity  which  is  our  goal.  To  obtain  inference  about  small  chemical  changes,  we  use  an  inverse 
problem  model  based  approach  on  the  data  that  we  collected. 
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Single  fiber  w /  coating  (10-8  m) 


FIGURE  1.  An  example  diagram  of  the  multiscale  geometry  of  ceramic  matrix  composites. 


40  nm 


FIGURE  2.  Diagram  of  the  degradation. 

EXPERIMENTAL  DESCRIPTION 

IR  spectroscopy  was  performed  on  silicon  nitride  carbon  based  CMC  samples  before  and  after  heat  treatment.  MIR 
measurements  were  made  using  a  Bruker  Vertex  80v  FTIR  attached  to  a  Harrick  Omni-Spec  probe.  The  probe  made 
surface  contact  where  45-degree  reflection  measurements  were  collected  and  compared  to  a  mirror  reference.  The 
samples  were  15.7  x  1.3  x  0.2  cm  in  dimension.  Each  sample  was  divided  along  its  length  into  1 1  blocks  of  equal  area. 
Prior  to  measuring  a  sample,  reference  intensity  data  was  obtained,  and  then  an  intensity  measurement  was  collected 
in  the  approximate  center  of  each  of  the  1 1  blocks.  This  procedure  was  considered  one  trial,  and  was  repeated  for  the 
second  and  third  trials.  The  reflectance  data  was  calculated  as  the  ratio  of  the  intensity  measurement  from  the  sample 
to  the  intensity  measurement  of  the  reference. 
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TABLE  1.  Acquisition  parameters. 

|  #  of  pts  |  1st  Pt  (1/cm)  |  2nd  Pt  (1/cm)  |  Resolution  |  Sample  scans 


Baseline 

Sample  13 

5807 

5998.352276 

399.247241 

8 

16 

Sample  16 

5807 

5998.352276 

399.247241 

8 

16 

Sample  4 

14787 

7499.386893 

369.834099 

4 

16 

Heat  10H 

Sample  16 

14787 

7499.386893 

369.834099 

4 

16 

Sample  32 

14787 

7499.386893 

369.834099 

4 

16 

Sample  1 

14787 

7499.386893 

369.834099 

4 

16 

Sample  13 

5807 

5998.352276 

399.247241 

8 

16 

Sample  29 

14787 

7499.386893 

369.834099 

4 

16 

Baseline  data  was  collected  on  samples  13  and  16  using  the  following  specifications:  5807  points  over  a  frequency 
range  of  5998.35  to  399.25  cm-1  with  a  resolution  of  8  cm-1. 

Three  samples  of  SiNC  (samples  4,  16,  and  32)  underwent  heat  treatments  of  1200  deg.  C  for  10  hours.  MIR  data 
was  collected  on  these  samples  using  the  same  method  as  before  except  under  the  following  specifications:  14787 
points  over  a  frequency  range  of  7499.39  to  369.83  cm-1  with  a  resolution  of  4  cm-1. 

In  addition,  three  samples  of  SiNC  (samples  1  and  29)  were  exposed  to  heat  treatments  of  1200  deg.  C  for  100 
hours.  These  samples  were  inspected  in  the  same  manner  as  samples  4,  16,  and  32.  Sample  13  was  exposed  to  the 
same  heat  treatment  as  samples  1  and  29  but  was  inspected  using  the  same  specifications  used  to  collect  the  baseline 
data.  A  summary  of  the  collected  data  is  given  in  Table  1 . 


DATA  ANALYSIS 


To  relate  changes  in  data  to  such  small  physical  changes,  we  need  to  analyze  the  data  with  models  that  give 
reflectance  from  parameters  that  relate  to  chemical  changes.  We  would  like  these  to  result  in  information  about  the 
underlying  amounts  of  chemicals  and  their  phase.  An  example  of  the  changes  in  the  data  is  given  in  Figure  3 

We  first  use  standard  optical  results  to  model  the  reflectance  and  the  refractive  index  of  the  material.  Throughout 
this  description,  we  use  the  notation  in  Table  2.  The  two  polarizations  of  the  reflectivity  are  given  by 


iia  cos  di  —  ( Pi2 ( x ; q^)  —  n\  sin2  0/) 

r±(x\q(N))  = - - - — 

iia  cos  0/  +  ( n 2 (x; gA) )  —  n\  sin2  0/)  2 

(1) 

r{x.q[N)] = ha  0 1  -  sin2  6/) 2  ~ n (*;  q(N)  ]  cos  6l 

nA  ( 1  -  R2{"Aqm) sin2 6l)  2  +  "(*; ?(iV) ) cos 

(2) 


which  is  derived  using  Snell’s  law  and  Fresnel’s  law  [8].  The  total  reflectance  is  given  by  the  the  average  of  the 
two  polarizations  of  the  reflectance.  These  are  the  moduli  squared  of  the  complex  two  polarizations  of  the  complex 
reflectivity: 

R{x\qW)  =  \  (|r±(x;^)|2+|r||(^;?W)|2)  ,  (3) 

where  |  •  |  is  the  modulus  of  the  complex  number  (z  =  x  +  iy  then  |  z|  =  x2  +y2)  [8] . 

The  refractive  index  of  a  material  depends  on  its  complex  dielectric  constant  through  the  following  relation: 

h2(x)  =  e2(x). 
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Reflectance  data 


Example  Data 


x  cm 

FIGURE  3.  Example  data  from  randomly  selected  data  sets. 


TABLE  2.  Notation. 

Symbol  |  Parameter  |  Value 


e, 

Incident  angle  of  beam 

45  degrees 

nA 

Refractive  index  of  surroundings 

1 

X 

Wave  number 

Dependent  variable,  varies  from  400  — 
1400  cm-1 

Dielectric  constant  or  permittivity  (complex  number) 

Varies  with  wave  number 

n 

Refractive  index  of  material  (complex  number) 

Depends  on  £ 

Perpendicular  polarization  reflectance 

Depends  on  h ,  Snell/Fresnel 

*11 

Parallel  polarization  reflectance 

Depends  on  n,  Snell/Fresnel 

* 

Total  Reflectance 

Average  of  R±  and  Sy 

We  model  the  dielectric  constant  as  a  function  which  is  based  on  assumptions  about  molecular  motions.  There  have 
been  many  models  proposed  for  the  dielectric  constant  which  vary  regarding  form,  and  number  of  variables.  We  will 
focus  on  models  that  use  springs  (or  oscillators)  to  model  molecular  vibrations.  We  will  denote  e(x)  the  complex 
dielectric  function  with  real  and  imaginary  parts  £1  and  £2,  respectively, 

£  (x)  =  £1  (x)  +  ie 2  (x) .  (4) 

The  dielectric  constant  has  dependent  variable  wave  number  x  cm-1.  The  classical  oscillator  model  is  given  by 


£(x) 


>+L 


Ml  —  X2)  —  iTj. 


J= 


(5) 


where  N  is  the  number  of  oscillators,  (Qj  is  the  resonant  frequency,  Fj  is  the  broadening,  and  Aj  is  the  product  of  the 
strength  and  the  number  of  molecules  of  that  oscillator  type.  This  is  known  as  the  Lorentz  model  in  dispersion  analysis 
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and  is  most  appropriate  for  crystalline  systems  where  the  environment  is  the  same  for  oscillators  of  the  same  kind  [9]. 
Efimov  proposed  a  model  where  the  resonant  frequency  is  normally  distributed  and  broadening  is  the  same  for  all  the 
oscillators  [1].  When  the  standard  deviation  is  much  larger  than  the  broadening,  this  model  can  be  approximated  by 


E2W  =  £ 


H 1  V27iAcOj 


-<°ojY 

I 


£1  (x)  =  £oo  +  H.  T.E2  (x) , 


(6) 


where  H.T.  refers  to  the  Hilbert  transform  [10,  9]. 

We  analyzed  our  data  sets  with  these  two  models.  For  N  6  {1,2,3, . .  .20},  we  estimated  the  parameters  in  the  above 
models  by  minimizing  the  residual  sum  of  squares.  We  used  Aikake  information  criterion  to  choose  between  Lorentz 
and  Gaussian  models.  For  almost  all  data  sets  and  oscillator  numbers,  the  curve  fit  was  better  for  the  Gauss  model  than 
the  Lorentz  model.  We  present  the  results  for  Gaussian  model  in  the  remainder  of  this  manuscript  [11]. 


RESULTS 

The  models  that  we  chose  do  not  appear  to  fit  the  data  especially  well.  Indications  of  this  were  observed  in  poor 
fits  of  the  curves,  many  local  minima  which  caused  the  minimization  of  the  residual  sum  of  squares  to  be  difficult  (we 
used  implicit  filtering  [12]),  and  the  large  variation  in  parameter  fits  across  (a)  samples,  (b)  locations  and  (c)  trials.  We 
have  presented  histograms  of  the  estimated  parameters  in  Figures  4-6  (a)-(c).  There  is  a  large  variation  of  the  estimated 
parameters  for  all  of  the  estimated  parameters.  Often  in  spectroscopy,  the  estimated  resonant  frequencies  ( (Boy’s)  are 
used  to  identify  chemical  species  and  the  amplitudes  (Ay’s)  are  used  to  determine  the  amount.  In  Figures  4  (a)-(c),  we 
see  that  there  are  regions  of  concentration  of  resonant  frequencies  but  exactly  how  many  there  are  is  not  discernible. 
There  does  appear  to  be  a  change  in  resonant  frequency  between  the  baseline  samples  (in  Figure  4(a)),  and  the  treated 
samples.  The  randomness  of  these  make  the  exact  differences  difficult  to  discern  and  impossible  to  relate  to  exact 
chemical  phases.  There  is  a  large  variation  in  the  estimates  of  the  parameters  A j  and  A C0j  though  there  are  less  regions 
of  high  concentrations  as  in  the  estimates  of  (Boy. 


FUTURE  WORK 

The  results  of  the  our  data  analysis  did  indicate  differences  between  the  different  treatment  levels.  The  apparent 
randomness  of  the  estimates  themselves  hindered  interpretation  of  these  changes.  There  could  be  several  causes  of 
this  issue.  We  could  be  fitting  model  discrepancy  with  extra  oscillators.  We  will  first  try  to  relax  the  assumption  that 
Tj  »  AcOj  that  we  made  in  going  from  the  Efimov  model  to  the  Gaussian  model  (6)  (as  in  [9]).  This  model  has  been 
validated  for  different  kinds  of  glass  which  leads  us  to  believe  it  will  fit  our  data  better  and  give  us  easier  to  interpret 
results. 
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FIGURE  4.  Histogram  of  too,-  for  100  hour  treated  data  set. 
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FIGURE  5.  Histogram  of  log10  Aj  for  100  hour  treated  data  set. 
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